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■ ^ 1 INTRODUCTION 



ABSTRACT 

We present wide field-of-view near-infrared imaging from the NTT and very deep op- 
tical imaging from the HST of the young merging galaxy NGC 3597. The morphology 
of the galaxy and the properties of the newly formed proto-globular clusters (PGCs) 
are examined. Our K band data reveals the presence of a second nucleus, which pro- 
vides further evidence that NGC 3597 is the result of a recent merger. Combining new 
K band photometry with optical photometry, we are able for the first time to derive 
a unique age for the newly formed PGCs of a few Myrs. This is consistent with the 
galaxy starburst age of < 10 Myrs. From deep HST imaging, we are able to probe the 
luminosity function ~ 8 magnitudes fainter than normal, old globular clusters, and 
confirm that the PGCs have a power-law distribution with a slope of ~ -2. 
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In the 1970's, the rich globular cluster (GC) populations 
seen around elliptical galaxies were used to argue against the 
idea that ellipticals formed from the simple merger of two 
spirals (van den Bergh 1975). The suggested solution to this 
problem was the creation of new GCs from the gas associated 
with the progenitor galaxies (Schweizer 1987; Ashman & 
Zepf 1992). The Hubble Space Telescope (HST), with its high 
spatial resolution, has indeed detected proto-globular clus- 
ter candidates in several merging galaxies. Over a dozen such 
systems have now been observed with HST (e.g. Holtzman 
et al. 1992; Whitmore et al. 1993; Whitmore & Schweizer 
1995; Schweizer et al. 1996; Holtzman et al. 1996; Miller et 
al. 1997; Zepf et al. 1999; Carlson et al. 1999). These stud- 
ies suggest that all mergers involving gaseous systems create 
proto-globular clusters (PGCs). However there are several 
outstanding issues concerning these PGCs, such as the total 
number created, their luminosity (mass) function, their de- 
struction rate and the overall specific frequency of the final 
system (e.g. Kissler-Patig, Forbes & Minniti 1998; Brodie 
et al. 1998; Zepf et al. 1999). 

Perhaps the first detailed study of PGCs in a merg- 
ing system was that of Lutz (1991). Lutz presented opti- 
cal imaging and spectroscopy of NGC 3597 (AM 1112-232). 



Classified as an SOpec, it reveals an extended structure with 
plumes but the main body of the galaxy resembles an r^' *- 
like surface brightness profile. Using the ESQ 2.2m telescope, 
under ~1.3 seeing conditions, Lutz detected and measured 
photometry for 31 unresolved objects around the galaxy. 
About half of these were suggested to be young, blue PGCs. 
Following this ground-breaking work by Lutz, Holtzman et 
al. (1996) re-observed NGC 3597 using the WFPC2 camera 
on HST. They detected over 70 PGC candidates. The bulk of 
these had similar colours suggesting a single age and metal- 
licity population. However from their V-R colours alone 
they were unable to identify a unique age for the PGCs. 

Here we present K band (2.2^m) imaging of NGC 3597 
taken under photometric conditions with ~0.6 seeing. The 
near-infrared has the advantage of being much less affected 
by dust extinction (Ak ~ 0.1 Ay). We also incorporate deep 
B and R band data from the HST archive. The near-infrared 
photometry, when combined with optical magnitudes, pro- 
vides a powerful constraint on the age of PGCs. After com- 
pleting our study, Carlson et al. (1999) published a study of 
NGC 3597 using the same deep HST imaging. However as 
described below, the deep B and R data alone do not pro- 
vide a significant improvement over the Holtzman et al. V 
and R band data in determining the age of the PGCs. Our 
main contribution is the addition of new K band photometry 
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which aUows the age of the PGCs to be uniquely constrained 
for the first time. The K band images also reveal a second 
nucleus, not obvious in the shorter wavelength HST images, 
but apparent from radio imaging. 



2 NEAR-INFRARED DATA 

2.1 Observations and Data Reduction 

K short (Ks; A/AA = 2.162/0.275) images of NGC 3597 were 
taken with the infrared camera SOFT on the ESO NTT over 
2 nights in 1999 March. The data were obtained under pho- 
tometric conditions and excellent seeing (0.6 ). SOFT has a 
field-of-view of ~ 5 x 5 and a pixel size of 0.292 . A total 
of 220 X 30sec exposures were taken to give a total effective 
exposure time of llOmin. Between each SOsec exposure, the 
telescope position was offset by ~1 . Standard stars from 
the list of Persson et al. (1999) were also taken before and 
after each block of galaxy observations. 

Flat-fields were made by median combining the stan- 
dard star images. These were determined to be superior to 
the dome and twilight flats. After subracting dark currents 
and flat -fielding, the individual images were carefully shifted 
and median combined to form a single final Ks image of NGC 
3597. The sky background of the final image is fiat to ^^2%. 

From our standard star images taken directly before 
and after the galaxy, we have determined a Ks zero point of 
22.23 ± 0.03, where the uncertainty represents the errror on 
the mean from several different images. We will not attempt 
to transform our Ks magnitudes to standard K band values, 
but simply note that they are similar and have been defined 
so that Ks = K for AO standard stars (Persson et al. 1999). 
Hereafter we simply use K to refer to our K short band data. 



2.2 Galaxy Modelling 

We have modelled the galaxy using the STSDAS task ellipse. 
The central position and the position angle of the model 
were fixed, but the ellipticity was allowed to vary. Stars and 
bright PGCs were masked out. We also used sigma clipping 
to exclude the most deviant data values from the model fit. 
Subtraction of the resulting galaxy model showed that it was 
a reasonable representation except in the very inner regions. 
The K surface brightness profile for the model is shown in 
Fig. 1. The outer regions are close to an r^'* surface bright- 
ness profile, as found in most elliptical galaxies. The total 
magnitude of the galaxy, from a curve-of-growth analysis, is 
K = 11.75 ± 0.05. We later use the model-subtracted image 
for the detection of globular clusters. 



3 RESULTS AND DISCUSSION 

3.1 The Host Galaxy 

The basic properties of NGC 3597 are summarised in Table 
1. Coordinates of the galaxy centre are taken from an HST 
R band (F702W) image, but are probably only accurate to 
within ±1 due to dust. However this position is consistent 
with the position of the brightest radio source (van Driel et 
al. 1991). The distance of 49 Mpc is based on a Virgocentric 
infiow corrected velocity of 3513 km s^^ and Ho = 75 km 



s~^ Mpc~^. At this distance 1 corresponds to 240 pc. The 
galaxy is fairly luminous at all wavelengths, and is undergo- 
ing a vigourous starburst (e.g. Lutz 1991; Kim et al. 1995; 
Rephaeli et al. 1995). 

In Fig. 2 we show the inner regions of our K image of 
NGC 3597. We also identify the location of the two radio 
sources mapped at 6cm by van Driel et al. (1991), which 
correspond spatially with the bright central source and a 
second source in our K image. The galaxy inner regions re- 
veal that the bright central source is connected via a 'bridge' 
to the second source about 3.8 (0.9 kpc) to the West. This 
could be simply a super starcluster, but it is not particuarly 
bright at optical wavelengths. At radio wavelengths, it ap- 
pears as a connected second source with a similar fiux (Seem 
= 18.0 mjy compared to the central source of 29.1 mjy). 
Thus it appears quite likely that the central source and the 
second source, seen in the K band and at 6cm, are the nu- 
clei of the merging galaxies. The dynamical timescale for 
these nuclei to merge is less than 10^ yrs, so we appear to 
be witnessing NGC 3597 in the very final stages of nuclear 
coascelence. 

With a projected nuclear separation of 0.9 kpc, NGC 
3597 is at a similar evolutionary stage to NGC 3256 (Norris 
& Forbes 1995) and NGC 6240 (van der Werf et al. 1993) 
with separations of 0.73 and 0.90 kpc respectively. The ra- 
tio of far-infrared to H2 mass is a measure of star formation 
efficiency, and has been shown to increase as the nuclear 
separation of two merging galaxies reduces to zero (Gao & 
Solomon 1998; Georgakakis, Forbes & Norris 1999). NGC 
3597 has a Ijfir/Mh^ of 25.6 (Wiklind et al. 1995), com- 
parable to 20 for NGC 3256 and 30 for NGC 6240. This 
further supports the case that NGC 3597 is near the nuclear 
coascelence stage of a merger sequence. 

Beyond our faintest isophote, deep optical images reveal 
plumes {^R ~ 25 mag arcsec"^) to the NW and SW (Lutz 
1991). As NGC 3597 is relatively isolated, the outer mor- 
phological disturbance is almost certaintly due to a merger 
of some sort. Indeed the morphology of the plumes resem- 
bles the merger simulation shown in Barnes (1998) at 300 
Myr after pericentre, and just before nuclear merger. 

The galaxy is currently undergoing a starburst. Optical 
spectra (e.g. Lutz 1991; Kim et al. 1995) reveal emission lines 
indicative of HII regions. However, as well as the current star 
formation activity, there is also evidence that the burst is 
not instantaneous but has proceeded for some time. Optical 
spectra clearly show weak Hydrogen absorption lines (along 
with the Hydrogen emission) indicative of an earlier phase 
of star formation. The radio spectral index is q = -0.84 
indicating that the emission is not dominated by HII regions 
but rather is due to non-thermal synchrotron emission from 
SNRs (Smith & Kassim 1993). Both the current and post 
starburst contribute to the far-infrared luminosity of L_f/h 
~ 7 X 10^" L0 which is consistent with the well known far- 
infrared vs radio correlation. Kim et al. (1995) estimated an 
H/? absorption EW of 2.5 A from their optical spectrum. 
Such an EW occurs in the first 10 Myrs of a starburst (e.g. 
Bruzual & Chariot 1993). This indicates that the starburst 
is very recent, starting less than 10 Myr ago and continuing 
to the present day. We compare this age estimate to that of 
the PGCs in the next section. 

Wikland et al. (1995) did not detect HI gas in NGC 
3597 but only placed a relatively high upper limit of yini < 
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2.6 X 10^° Mq. Based on the FIR luminosity and a Galactic 
gas-to-dust ratio, Lutz (1991) estimated an HI gas mass of 
3 X 10* M0, i.e. almost a factor of one hundred below the 
current Wikland et al. limit. It would be interesting to place 
tighter limits on the HI mass and investigate the possibility 
that atomic hydrogen gas is being compressed during the 
merger into molecular gas. 

Mihos & Hernquist (1994) have modelled the star for- 
mation history of equal mass spiral galaxies. They found 
that for systems without strong bulges (e.g. Sc spirals) , most 
of the merger-induced star formation occurs at the pericen- 
tre encounter over a timescale of a few Myrs. However the 
presence of a bulge inhibits gas flow and delays the main 
starburst until nuclear coalescence. This starburst also lasts 
only a few Myrs. The age of the starburst in NGC 3597 (i.e. 
< 10 Myrs) and the galaxy morphology (i.e. two nuclei and 
a common r^' * like envelope) suggests that the latter situ- 
ation is the correct one. Thus it appears that NGC 3597 is 
likely the product of two near equal mass spirals, at least 
one of which contained a strong bulge. Furthermore we are 
currently witnessing it at the nuclear merger stage, which is 
also the period of its most intense star formation activity. 



3.2 The Proto-Globular Clusters 

Our K image of NGC 3597 reveals a number of bright un- 
resolved sources. The vast bulk of these are PGCs. They 
were first noticed by Lutz (1991) and later re-observed with 
HST (Hohzman et al. 1996; Carlson et al. 1999). In the dis- 
cussion that follows we will assume a Galactic extinction 
towards NGC 3597 of Ay = 0.12 (As = 0.16, Ar = 0.08, 
and Ak = 0.01) as used by Holtzman et al. Although the 
galajcy shows a dust lane to the North of the nucleus, the 
internal extinction towards most of the PGCs is probably 
not large given the relatively uniform V-R colors noted by 
Holtzman et al. (see their figure 6). 

Lutz (1991) tabulated 31 sources. From Holtzman et al. 
(1996), it appears that 9 of them are bona fide PGCs, i.e. 
numbers 10, 11, 13, 14, 17, 19, 21, 23 and 24 from his table 
(the other 22 sources are either foreground stars or back- 
ground galaxies). The mean B-V colour for these 9 sources, 
corrected for Galactic extinction, is 0.40 with a la error on 
the mean of 0.03. This colour and range for the PGCs is 
shown in Fig. 3, along with the evolutionary track for a so- 
lar metallicity, Salpeter IMF single stellar population from 
Bruzual & Chariot (1993). The IMF and metallicity of the 
PGCs are of course unknown, so any derived ages will be a 
guideline only. The average PGC colour intersects the track 
at about 6-10, 100 and 1000 Myrs. Thus it is impossible to 
derive a unique age from the B-V colour alone. 

Hohzman et al. (1996) used WFPC2 to obtain F555W 
(V) and F702W (R) photometry of the PGCs. They detected 
72 PGCs in the PC chip down to V ~ 25.5. The first source 
listed by Holtzman et al. (Holtz ID = 1) is probably a star 
with V = 20.04 and V-K = 1.4 which suggests it is a GO 
dwarf star. For the 60 globular clusters with V < 25, the 
mean colour and error on the mean is V-R = 0.34 ± 0.03 
(again corrected for Galactic extinction). This is shown in 
Fig. 4 along with the Bruzual & Chariot (1993) evolutionary 
track. As with B-V, the V-R colours do not give a single 
age but several possibilities, i.e. ~ 6, 12, 80, 200 Myrs. The 



V-R colours do however rule out the oldest B-V age of 1 
Gyr. 

Using the HST PGC positions, we measured K magni- 
tudes using daophot. This was carried out on the galaxy sub- 
tracted image, using small aperture magnitudes for which we 
then applied an aperture correction (based on objects away 
from the galaxy centre). We detected 31 PGCs in our im- 
age down to K ~ 22.8, although with large errors at faint 
magnitudes. Our K magnitudes and measurement errors are 
listed in Table 2, along with Holtzman et al. ID number and 
offset from the galaxy centre. All of these PGCs lie within 
~ 20 of the galaxy centre (i.e. detected in the PC chip by 
Holtzman et al. ). 

In Table 3 we list the V magnitude, V-R colours and 
errors from Holtzman et al. (1996). The errors are photo- 
metric measurement errors, and do not include zeropoint 
errors (on the order of 0.05-0.1 mag.). These errors have 
been combined in quadrature to obtain the final V-K error. 
The mean colour and error on the mean is V-K = 0.95 ± 
0.14. This range of colours is shown in Fig. 5, along with 
the Bruzual & Chariot (1993) evolutionary track. The extra 
'leverage' from the K band means that the V-K colours of 
the PGCs only intersect the track at one place, i.e. around 
5 Myrs. This is consistent with the youngest age suggested 
by the B-V and V-R colours, and rules out the older ages 
at the 3a" level. Any correction for internal extinction, or if 
the PGCs had supersolar metallicities, would make the age 
even less than 5 Myrs. Such an age should be regarded as 
somewhat qualitative since it is depedent on unknown fac- 
tors (e.g. the IMF) and on the particular stellar population 
model used (i.e. Bruzual & Chariot 1993). 

We have obtained deep F450W (B) and F702W (R) 
HST images of NGC 3597 from the CADC HST archive. 
The 8 images total 5200sec in B and 5000sec in R, and cover 
essentially the same area as the shallower Holtzman et al. 
(1996) data. They were average combined using the STS- 
DAS task gcombine, which effectively removes cosmic rays. 
The high throughput of the F702W filter and the spectral 
energy distribution of globular clusters means that the R im- 
age will be somewhat deeper than the B image. Nevertheless 
both images should reach magnitudes typical of the brighter 
GCs associated with the original, old population from the 
progenitor galaxies. The previous HST imaging published 
by Holtzman et al. (1996) had an F702W exposure time of 
llOOsec and was thus not deep enough to reach the old pop- 
ulation and only probe the bright end of the new population. 
This data was also recently used by Carlson et al. (1999), 
and reaches about 1.6 mags deeper in R than the Holtzman 
et al. study. 

In order to locate sources, we used daofind with a con- 
servative 4(T per pixel detection criterion. For the PC chip, 
we fit the galaxy isophotes using ellipse and subtracted off a 
model before running daofind. Because of the confused na- 
ture in the central regions in the B image, due to dust and 
young star formation, we have excluded the central 9 radius 
from the automatic detection. Photometry on all four chips 
was measured using phot with a 2 pixel radius aperture. 
Aperture corrections to 0.5 and zero points were taken from 
Holtzman et al. (1995, 1996). The objects on the PC chip 
are marginally resolved. Holtzman et al. (1996) showed that 
the bulk of objects required a typical correction of 0.22™, 
in addition to the above aperture correction, to include all 
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of the light. Variations with the size of the object and its 
position on the chip, could introduce another source of error 
on the order of 0.1™. However the correction to total light 
affects the two filters almost equally, so there is little error 
introduced in the final B-R colour. 

The resulting object lists were then checked visually on 
the image display to exclude the very few remaining ob- 
vious bright stars and background galaxies. At this stage, 
our candidate GC list contained 292 sources from all four 
chips (excluding the central 9 of the PC chip). A colour- 
magnitude diagram is shown in Fig. 6. We have adopted 
selection cuts of -0.5 < B-R < 2.0 and B < 27. The colour 
cuts correspond to the full range expected of GCs older than 
a Myr (Bruzual & Chariot 1993) . The faint magnitude limit 
was chosen to avoid any colour bias in our sample. Within 
these selection criteria we have 239 objects, the vast bulk of 
which will be bona fide GCs (as our contamination rate is ^ 
5%). 

A histogram of B-R colour is shown in Fig. 7. The GCs 
have an average colour of B-R = 0.66 ± 0.03 (error on the 
mean). The spread in colour is consistent with photometric 
errors suggesting that the PGCs are close to a single age 
and metallicity population. This mean value is shown on the 
evolutionary tracks of Bruzual & Chariot (1993) in Fig. 8. 
From B-R colour alone the GCs may have several possible 
ages, but from our K band imaging above it is clear that 
the correct age is ~ 5 Myrs. Thus the majority of detected 
objects in our deep HST images are PGCs. 

Have we detected any GCs from the progenitor galaxies 
? Milky Way GCs have a mean B band luminosity function 
that is roughly Gaussian with a peak at Ms ~ -6.6, with the 
brightest GC {uj Cen) having Mb = -10.7. These correspond 
to B = 26.85 and 22.75 respectively, at the distance of NGC 
3597. Uncertainty in the distance makes these magnitudes 
uncertain by about ± 0.2™. Milky Way GCs have an average 
extinction-corrected colour of (B-R)o ~ 1.2. Examination 
of the colour-magnitude diagram reveals 32 GCs in the ex- 
pected colour and magnitude range (shaded region in Fig. 
6) , although perhaps half of these will simply be PGCs with 
apparently red B-R colours due to photometric errors. If we 
assume that say 20 are true old GCs then we can crudely 
estimate the total population of progenitor GCs. The HST 
images cover about half of the area out to a galactocentric 
radius of 125 arcsec, suggesting 40 GCs within this radius. 
Our HST data only reach to B ~ 26.5, which is slightly less 
than the expected peak at B = 26.85. This suggests that 
we are sampling about 40% of the GC population in mag- 
nitude terms. Thus another correction of ~2.5x is required, 
giving a total old GC population of about 100. This crude 
calculation could easily be in error by a factor of two, but is 
not vastly different to the halo GC population of the Milky 
Way (i.e. ~120). Unfortunately such a calculation does not 
constrain the progenitor types. Deeper B band observations, 
over a somewhat wider field of view, would provide better 
constraints on the original GC systems. 

Globular cluster luminosity functions in old ellipticals 
show a Gaussian log-normal distribution, with a peak, 
or turnover, magnitude and fewer GCs at low luminosity 
(mass). However, for the PGC systems studied to date the 
luminosity function does not peak but continues as a power- 
law down to low luminosities. Published slopes include -1.78 
(NGC 4038/9; Whitmore & Schweizer 1995), -2.1 (NGC 



3921; Schweizer et al. 1996), -1.84 (NGC 7252; Miller et 
al. 1997) and -1.8 (Zepf et al. 1999). For the closer systems, 
the PGCs in the PC chip are marginally resolved giving 
sizes consistent with GCs rather than open clusters (which 
also have a power-law distribution). Recently Carlson et al. 
(1999) have examined the luminosity function for NGC 3597. 
They found a slope of -2.0 to be a reasonable representa- 
tion. Here using the same HST data, we also probe the R 
band luminosity function. 

Before probing the GC luminosity we need to esti- 
mate the incompleteness of our R band detections. This was 
achieved by simulating GCs using the addstartask, and mea- 
suring the detected fraction as a function of R band magni- 
tude. An actual WFC chip was used to reproduce the correct 
noise characteristics. Care was taken to avoid any blending 
of artifical GCs. The resulting completeness function for the 
WFC chips is shown in Fig. 9. Our 50% incompleteness level 
is at R ~ 26.8. In Fig. 10 we show the R band luminosity 
function in log-log space, for the three WFC chips and cor- 
rected for incompleteness. It is very similar to that shown 
in figure 10 of Carlson et al. (1999). A Milky Way like GC 
system with a peak magnitude of M_r — -8, would have an 
apparent magnitude of R ~ 25.5 at the distance of NGC 
3597 and be an additional 5.5 magnitudes brighter if they 
were 5 Myr old instead of 12 Gyr old (assuming a solar 
metallicity population). This gives an expected peak of R 
~ 20, with the faint (lower mass) limit of R ^ 23. By con- 
strast, the observed PGC luminosity function continues to 
rise to faint magnitudes (R --^ 28) with a power-law like 
slope of around -2. At the faint end of our luminosity func- 
tion, assuming a Salpeter IMF, derived masses correspond 
to ~ 200 Mq. At 8 magnitudes fainter than the expected 
peak, these data probe deeper than any previously published 
study of PGCs. We note that similar mass calculations by 
Carlson et al. (1999) assumed an age of ~500 Myrs, which 
means that their masses are overestimates. Thus, like other 
merging galaxies, the low mass PGCs in NGC 3597 must 
be destroyed (e.g. by tidal disruption or evaporation) over 
time if their luminosity functions are to resemble those of 
old ellipticals (see Zepf et al. 1999). 

Elmegreen & Efremov (1997) have proposed that a uni- 
versal mechanism exists for GC formation. In their model, 
GCs form in high pressure regions with an assumed power- 
law mass distribution of slope '^ -2. They suggest that a 
GC system is formed initially without a characteristic mass, 
but as low mass GCs are destroyed over time a characteris- 
tic mass develops. Over ~ 10 Gyr, the luminosity function 
grows to resemble that of the Milky Way GC system with a 
characteristic mass corresponding to M_r ~ -8.0 

In our K image there are many unresolved sources be- 
yond the central ~ 20 (i.e. the region covered by Holtzman 
et al. 1996). We used daofind to find all sources 3a above 
the background noise. After rejecting 3 obvious bright stars, 
1 galaxy and objects fainter than with K = 23, we were left 
with a list of 142 candidate PGCs. These objects are listed 
in Table 4 along with their K magnitudes and photometric 
errors. Some of the brighter objects will be foreground stars, 
and some will be compact background galaxies. Indeed, after 
examining the area in common with the deep HST images, 
we have excluded a further three stars and two galaxies (as 
noted in Table 4). As the HST images do not cover the whole 
area of our K image, we are unable to confidently remove 



© 0000 RAS, MNRAS 000, |-?? 



NGC 3597 Globular Clusters 



all stars and galaxies. In the absence of overlapping HST 
data (or optical magnitudes), it is difficult to make conclu- 
sive statements about these outer objects, but many will be 
PGCs. 



4 CONCLUDING REMARKS 

The combination of near-IR and optical imaging from the 
NTT and HST telescopes has provided new insights into the 
merging system NGC 3597. For the host galaxy we confirm 
an r^'* like surface brightness profile in the outer regions 
and discover the presence of two closely separated (Ar ~ 
0.9 kpc) nuclei. The two nuclei seen in our K band image 
correspond to those seen at radio wavelengths and are prob- 
ably the nuclei of the progenitor galaxies at the final stages 
of coalescence. Various properties of the galaxy starburst 
suggest that the burst occured less than 10 Myrs ago and 
may have involved at least one early type spiral. 

A number of proto-globular clusters (PGCs) are identi- 
fied in NGC 3597 with average colours of V-K = 0.95. The 
extra leverage provided by our K band photometry allows 
us to derive a unique age for the PGCs of a few Myr. From 
deep B and R HST imaging we detect ~ 300 globular clus- 
ters, the vast bulk of which are PGCs. We derive an R band 
luminosity function which reaches 8 magnitudes fainter than 
the expected characteristic mass of a standard old globular 
cluster luminosity function, and confirm the slope of the lu- 
minosity function to be about -2. A small number of glob- 
ular clusters have properties consistent with Galactic ones, 
and we speculate that they were associated with the pro- 
genitor galaxies. 
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Figure 1. K band surface brightness profile for NGC 3597. Errors are on the order of the symbol size. The surface brightness in the 
outer parts is close to an r^/** profile, i.e. a straight line. 
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Figure 3. B— V colour evolution for a single stellar population of solar metallicity from Bruzual &: Chariot (1993). The horizontal line 
and dashed lines show the mean B— V colour and la range for 9 proto-globular clusters (PGCs) from Lutz (1991). The PGCs could have 
several possible ages from ^ 6 Myrs to 1 Gyr. 
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Figure 4. V-R colour evolution for a single stellar population of solar motallicity from Bruzual & Chariot (1993). The horizontal line 
and dashed lines show the mean V-R colour and la range for 71 proto-globular clusters (PGCs) from Holtzman et al. (1996). The PGCs 
could have several possible ages from ~ 6 Myrs to 200 Myr. 
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Figure 5. V-K colour evolution for a single stellar population of solar metallicity from Bruzual & Chariot (1993). The horizontal line 
and dashed lines show the mean V— K colour and la range for 31 proto-globular clusters (PGCs) from this paper. The PGCs have an 
age of ~ 5 Myr. 
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Figure 6. Colour-magnitude diagram for the detected objects from deep WFPC2 data. The dotted line shows the region B < 27, 
—0.5 < B-R < 2.0 used to select globular clusters. Typical photometric error bars are shown on the left. The shaded region represents 
the apparent magnitude and colour range for Milky Way globular clusters placed at the distance of NGC 3597. The data show that 
the majority of globular clusters have B-R ~ 0.66, which are identified as proto-globular clusters. There are a few old. Milky Way like 
globular clusters. 
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Figure 7. B-R colour histogram of tlie 239 selected globular clusters from deep WFPC2 data. The majority of globular clusters have 
B-R f^ 0.66, which are identified as proto— globular clusters. 
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Figure 8. B— R colour evolution for a single stellar population of solar metallicity from Bruzual & Chariot (1993). The horizontal line 
and dashed lines show the mean B— R colour and Icr range for ^^239 proto-globular clusters (PGCs) from this paper. The age of the 
PGCs is consistent with ~ 5 Myr determined from V— K colours. 
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Figure 9. Completeness function for globular clusters from simulations. Brighter than R = 26, essentially all globular clusters are 
detected. 
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Figure 10. Logarithmic globular cluster luminosity function in the R band. The solid line shows the raw data, and the dashed line 
after correction for incompleteness. The objects brighter than R = 24 are not shown. The straight line shows a power-law slope of —2. 
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TABLE 1 
NGC 3597 Properties 



Proj^erty 



a (J2000) 11 14 42.0 1 

S (J2000) -23 43 40.0 1 

Type SOpec RC3 

B", 13.34'" 2 

B-V 0.67"" 3 

Distance 49.0 Mpc 2 

Mb -20.12"- 1 

MjT -21.71"" 1 

L.B 1.72 X 10"Lo 1 

Ljr 9.82 X 10 = Lo 1 

Ljrj.R 7.23 X 10i°Lo 2 

Ljc 1.2 X 10"= erg/s 4 

Laoc™ 2.7 X 10=8 erg/s S 

Mjy^ 2.82 X IO^Mq 2 

Mj^r < 2.6 X 10i°Mo 2 

M<,„„ 4.85 X lOi^Mo 2 

NOTES: 1 = this paper; 2 = Wiklind et al. (1995); 3 = Lutz 
(1991); 4 = Repliaeli et oZ. (1995); 5 = Smitli e< al. (1993) 
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TABLE 2 

Globular Cllister Ks band Photometry 



Holtz ID 






Ks 

(mag) 


<t(Ks) 
(mag) 


2 


-8.9 


-9.8 


19.86 


0.15 


3 


-11.2 


7.6 


20.21 


0.22 


4 


-3.6 


-17.0 


20.33 


0.21 


a 


11.2 


-0.7 


20.43 


0.16 


6 


-10-3 


0.3 


20.92 


0.47 


7 


8.3 


4.7 


19.94 


0.13 


8 


3.2 


-5.9 


20.84 


0.18 


9 


-4.7 


15.0 


20.11 


0.15 


10 


6.2 


-6.8 


21.12 


0.36 


11 


-7.0 


1.6 


21.05 


0.16 


12 


-7.4 


-2.2 


21.98 


0.33 


13 


-3.7 


4.4 


20.97 


0.26 


14 


16.9 


11.5 


21.55 


0.28 


15 


6.0 


-15.0 


21.61 


0.44 


16 


-5.7 


-1.7 






17 


-12.1 


-8.3 


21.78 


0.35 


18 


-12.6 


2.0 


21.95 


0.53 


19 


-3.2 


6.9 


21.28 


0.26 


20 


-5.0 


-5.2 


19.51 


0.13 


21 


8.1 


3.8 


20.02 


0.17 


22 


7.3 


-6.1 






23 


-6.1 


-0.6 






24 


-2.9 


12.2 






25 


-5.3 


-7.1 


21.70 


0.20 


26 


1.5 


-12.1 


22.38 


0.80 


27 


-3.3 


6.6 


21.39 


0.25 


28 


-10. 


-0.5 


22.72 


0.53 


29 


6.1 


7.7 


20.84 


0.31 


30 


0.3 


15.2 






31 


11.9 


-3.4 


22.11 


0.57 


32 


-0.2 


-14.7 






33 


4.8 


5.5 


22.40 


0.31 


34 


4.5 


-7.4 


22.31 


0.36 


35 


-0.4 


4.7 


21.82 


0.28 


36 


-11.4 


5.0 


22.76 


0.80 


37 


16.5 


9.9 






38 


5.5 


8.1 


22.13 


0.28 


39 


12.3 


5.5 


22.24 


0.62 


40 


-4.6 


-10.7 







NOTES: Proto— globular clusters v^ritliin ■ — ^20 of the galaxy cen- 
tre. ID number and offset from the galaxy centre is from Holtzman 
Bt al. (1996). Ks magnitudes and measurement errors are from 

this paper. 
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TABLE 3 

Globular Cllister Colours 



Holtz ID Lut 


z ID V 


o-(V) 


V-R 


<t(V-R) 


V-K 


<t(V-K) 


2 


13 20.31 


0.01 


0.30 


0.01 


0.46 


O.IS 


3 


21 20.39 


0.01 


0.21 


0.01 


0.18 


0.22 


4 


10 20.58 


0.01 


0.21 


0.01 


0.25 


0.21 


S 


17 21.08 


0.01 


0.31 


0.01 


0.65 


0.16 


6 


20.84 


0.01 


0.21 


0.01 


-0.08 


0.48 


7 


19 21.39 


0.02 


0.49 


0.04 


1.45 


0.13 


8 


21.24 


0.01 


0.27 


0.01 


0.41 


0.18 


9 


24 21.02 


0.01 


0.21 


0.02 


0.91 


0.15 


10 


21.62 


0.01 


0.26 


0.01 


0.50 


0.36 


11 


21.92 


0.02 


0.31 


0.04 


0.88 


0.16 


12 


21.88 


0.02 


0.29 


0.04 


-0.10 


0.32 


13 


21.66 


0.03 


0.29 


0.05 


0.69 


0.26 


14 


23 22.32 


0.02 


0.30 


0.02 


0.77 


0.28 


IB 


11 22.16 


0.02 


0.26 


0.01 


0.56 


0.44 


16 


22.32 


0.04 


0.57 


0.07 






17 


14 21.91 


0.02 


0.18 


0.03 


0.14 


0.35 


18 


22.14 


0.02 


0.26 


0.03 


0.19 


3.61 


19 


22.33 


0.02 


0.30 


0.03 


1.05 


0.26 


20 


22.00 


0.03 


0.91 


0.04 


2.48 


0.13 


21 


22.75 


0.18 


0.90 


0.21 


2.73 


0.25 


22 


22.98 


0.03 


0.42 


0.04 






23 


22.48 


0.06 


0.53 


0.08 






24 


22.68 


0.03 


0.23 


0.04 






26 


23.31 


0.03 


0.S3 


0.04 


1.61 


0.20 


26 


22.82 


0.03 


0.16 


0.04 


0.44 


0.80 


27 


22.92 


0.03 


0.28 


0.04 


1.53 


0.25 


28 


23.50 


0.15 


0.28 


0.22 


0.78 


0.55 


29 


23.50 


0.05 


0.19 


0.05 


2.66 


0.31 


30 


22.94 


0.04 


0.12 


0.06 






31 


22.99 


0.05 


0.15 


0.06 


0.87 


0.57 


32 


23.11 


0.04 


0.17 


0.06 






33 


23.11 


0.05 


0.27 


0.07 


0.71 


0.31 


34 


23.19 


O.OS 


0.20 


0.07 


0.87 


0.36 


3S 


23.83 


0.07 


0.92 


0.08 


2.00 


0.29 


36 


23.31 


0.05 


0.12 


0.07 


0.55 


0.80 


37 


23.68 


0.05 


0.30 


0.06 






38 


24.00 


0.11 


0.26 


0.33 


1.87 


0.29 


39 


23.71 


0.06 


0.14 


0.08 


1.47 


0.11 


40 


23.46 


0.05 


0.30 


0.06 
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TABLE 4 
Outer. Object Ks band Photometry 



Aq 


A<5 


Ks 


<t(Ks) 


Aq 


A,5 


Ks 


<t(Ks) 


Aor 


A<S 


Ks 


^(Ks) 


(") 


(") 


(mag) 


(mag) 


(") 


(") 


(mag) 


(mag) 


(") 


(") 


(mag) 


(mag) 


35.8 


-103.9 


20.59 


0.25 


60.5 


-38.2 


20.90 


0.20 


-3.8 


40.4 


21.45 


0.35 


-102.4 


-98.7 


21.10 


0.29 


31.6 


-37.7 


21.26 


0.27 


-10.2 


40.5 


20.64S 


0.18 


17.9 


-97.7 


21.41 


0.24 


77.2 


-36.9 


21.31 


0.27 


-101.9 


40.4 


21.97 


0.60 


-7.5.4 


-96.3 


19.83 


0.11 


82.1 


-36.6 


21.20 


0.4.5 


25. S 


41.9 


20.63 


0.23 


-8.5.2 


-96.1 


21.59 


0.42 


7.9 


-34.9 


22.54 


0.39 


73.8 


42.4 


22.81 


0.81 


-9.8 


-94.6 


21.21 


0.33 


-6.5 


-34.9 


22.65 


0.95 


-94.6 


42.6 


19.28 


0.13 


.37.0 


-94.0 


20.08 


0.16 


52.4 


-34.2 


22.30 


0.42 


86.9 


42.8 


22.59 


2.40 


97.9 


-91.9 


22.61 


0.42 


70.8 


-33.1 


18.63 


0.09 


28.3 


43.4 


21.69 


0.42 


89.8 


-90.6 


21.24 


0.31 


77.0 


-32.8 


22.89 


0.90 


-55.6 


44.7 


21.39 


0.34 


-30.6 


-88.6 


19.84 


0.12 


90.0 


-27.3 


19.24 


0.12 


-109.4 


45.1 


20.72 


0.23 


-68.0 


-88.2 


21.87 


0.79 


52.0 


-26.2 


18. .53 


0.09 


4.9 


45.8 


20.20 


0.18 


-46.2 


-86.7 


19.61 


0.14 


6.6 


-26.2 


21.48 


0.33 


90.3 


48.6 


19.06 


0.11 


-24.0 


-86.5 


21.45 


0.39 


73.0 


-22.6 


20.98 


0.27 


-18.0 


48.5 


19.64 


0.15 


12.9 


-85.7 


16.21 


0.03 


-59.3 


-18.6 


18.29* 


0.08 


96.2 


50.5 


19.60 


0.13 


104.1 


-82.9 


20.46 


0.16 


54.3 


-18.2 


20.43 


0.17 


-99.5 


50.9 


21.32 


0.26 


54.3 


-81.6 


21.46 


0.32 


-36.7 


-17.9 


22.02 


0.44 


-62.5 


53.2 


20.75 


0.2S 


58.8 


-81.4 


18.97 


0.10 


86.5 


-17.6 


21.01 


0.28 


56.3 


58.1 


19.84 


0.18 


34.0 


-81.2 


16.94 


0.04 


87.7 


-14.8 


21.36 


0.29 


-55.7 


59.3 


21.29 


0.24 


-21.5 


-81.2 


22.37 


0.33 


-51.2 


-12.4 


22.77 


0.45 


49.2 


62.3 


20.35 


0.19 


-104.2 


-80.5 


19.14 


0.11 


102.1 


-12.0 


19.14 


0.09 


-97.7 


62.5 


16.45 


0.03 


-21.1 


-78.8 


18.98 


0.11 


-96.0 


-8.6 


21.03 


0.25 


66.5 


66.9 


20.49 


0.22 


-103.1 


-78.4 


20.90 


0.21 


55.2 


-3.4 


20.13 


0.20 


20.0 


67.6 


22.69 


0.84 


37.4 


-78.2 


20.12 


0.17 


85.2 


-2.8 


19.98 


0.14 


-24.4 


69.0 


21.13 


0.24 


102.2 


-76.6 


22.66 


1.86 


-101.0 


-1.1 


21.11 


0.25 


-56.8 


69.9 


21.24 


0.29 


61.7 


-73.3 


21.24 


0.29 










-102.8 


72.1 


20.08 


0.14 


-2.9 


-71.2 


22.89 


2.87 


— - 


— - 


— - 


— - 


62.5 


74.3 


22.45 


0.75 


-38.0 


-70.7 


22.16 


0.36 










-32.3 


76.9 


21.41 


0.40 


28.3 


-70.6 


21.65 


0.46 


91.4 


3.3 


22.59 


0.63 


-38.8 


77.4 


21.98 


0.40 


-20.5 


-69.5 


22.38 


0.62 


-32.3 


11.2 


21.41 


0.43 


93.0 


77.8 


22.39 


1.12 


-91.8 


-67.6 


21.50 


0.43 


-42.0 


11.9 


21.19 


0.22 


35.4 


79.0 


20.95 


0.23 


44.1 


-67.2 


22.39 


0.48 


-91.9 


11.4 


18.90* 


0.10 


14.4 


80.9 


20.83 


0.18 


-60.8 


-65.4 


21.27 


0.28 


-41.6 


12.4 


20.76 


0.20 


1.9 


85.4 


21.18 


0.33 


93.0 


-64.3 


21.17 


0.30 


-34.6 


14.7 


21.91 


0.34 


-7.5 


90.4 


19.20 


0.12 


80.3 


-61.8 


19.19 


0.10 


61.9 


16.2 


19.95 


0.16 


100.5 


91.3 


21.15 


0.26 


45.1 


-61.1 


18.84 


0.10 


97.1 


17.2 


20.08 


0.12 


-3.7 


92.4 


19.53 


0.13 


-4.9 


-59. S 


22.70 


0.72 


8.5.2 


17.5 


22.01 


0.59 


-8.5 


92.6 


20.66 


0.23 


32.7 


-58.4 


18.77 


0.10 


32.9 


19.8 


19.70 


0.14 


-38.5 


94.4 


19.14 


0.12 


-82.5 


-56.4 


20.73 


0.25 


-55.0 


21.1 


21.26 


0.26 


-12.6 


94.5 


20.91 


0.28 


-91.3 


-55.6 


19.07 


0.11 


-36.6 


23.3 


20.70 


0.22 


-32.8 


95.0 


18.32 


0.08 


12.6 


-53.5 


20.17 


0.14 


65.2 


23.9 


20.96 


0.28 


-90.0 


95.4 


20.92 


0.52 


-70.6 


-53.6 


21.38 


0.40 


-45.5 


25.5 


21.24 


0.32 


52.8 


96.1 


20.01 


0.16 


25.2 


-52.0 


21.40 


0.39 


-55.2 


25.5 


18.59* 


0.09 


84.9 


98.7 


20.66 


0.24 


84.2 


-48.2 


21.52 


0.31 


-43.3 


28.2 


19.80 


0.16 


6.2 


104.0 


21.89 


0.39 


-4.2 


-47.6 


19.56 


0.12 


-64.6 


29.0 


22.18 


2.21 


-89.3 


105.3 


20.41 


0.18 


-93.9 


-47.2 


21.68 


0.34 


11.5 


31.1 


21.61 


0.45 


-96.4 


105.6 


17.52 


0.05 


-22.1 


-46.6 


19.24 


0.13 


21.1 


33.5 


21.67 


0.36 










-75.7 


-43.4 


20.99 


0.28 


-6.4 


34.3 


22.85 


0.85 










63.2 


-43.0 


22.81 


1.65 


63.8 


35.1 


19.80 


0.15 










-25.3 


-39.8 


21.15 


0.24 


-10.3 


36.9 


20.50B 


0.20 










37.8 


-38.5 


21.43 


0.39 


-6.3 


38.3 


21.79 


0.45 











NOTES: Proto— globular cluster candidates in Ks iniage beyond 20 galactocentric distance. RA and Dec. 
offset are from the galaxy centre. Objects ixiarked ^ are identified as galaxies and * as foreground stars from 
HST images. See text for details. 
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